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Abstract. Stability properties of magnetic-field configurations containing the toroidal and axial field are considered. 
The stability is treated by making use of linear analysis. It is shown that the conditions required for the onset 
of instability are essentially different from those given by the necessary condition d(sB v )/ds > 0, where s is the 
cylindrical radius. The growth rate of instability is calculated for a wide range of the parameters. We argue that 
the instability can operate in two different regimes depending on the strength of the axial field and the profile of 
the toroidal field. 
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1. Introduction 

Turbulence generated by MHD instabilities can play an 
important role in enhancing transport processes in various 
astrophysical bodies, such as accretion and protoplanetary 
disks, galaxies, stellar radiative zones, etc. The anomalous 
turbulent transport can be particularly important in mag- 
netized gas where a wide variety of MHD instabilities can 
occur (see, e.g., Barnes et al. 1999). In this case, the onset 
of instability can be caused both by hydrodynamic mo- 
tions (for instance, differential rotation; see, e.g., Velikhov 
1959; Chandrasekhar 1960) or unstable magnetic configu- 
rations. Which field strength and topology can sustain a 
stable magnetic configuration is still rather uncertain de- 
spite extensive work (see Borra et al. 1982; Mestel 1999 
for review). 

Most likely, the best-studied magnetic configuration is 
one with a purely toroidal field. Ever since the paper by 
Tayler (1973), it has been known that toroidal fields can 
be unstable close to the axis of symmetry, if there is a 
non-zero electric current density on the axis. The growth 
rate of this instability is expected to be of the order of the 
time taken for an Alfven wave to travel around the star 
on a toroidal field line. However, even a purely toroidal 
field is stable if it decreases rapidly with the cylindrical 
radius s. For instance, Tayler (1973; see also Chanmugam 
1979) argued that the toroidal field B v is stable against 
axisymmetric perturbations if it satisfies the condition 
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d(Btp/ s)/ds < and to non-axisymmetric perturbations if 
d(sB 2 )/ds < 0. Note that a purely toroidal field can also 
be subject to the magnetic buoyancy instability (Parker 
1955; Gilman 1970; Acheson 1978) but the Tayler insta- 
bility likely appears first as the strength of the toroidal 
field increases (Spruit 1999). 

The stability of a purely toroidal field in the radiative 
zones of stars and accretion disks has been studied by a 
number of authors. Numerical modeling by Braithwaitc 
(2006) confirms that the toroidal field with B v oc s or 
cx s 2 is unstable to the m = 1 mode (m is the azimuthal 
wave number) as predicted by Tayler (1973). The linear 
stability of the toroidal field in rotating stellar interiors has 
been considered by Kitchatinov & Riidiger (2007), who 
conclude that the magnetic instability is essentially three- 
dimensional and that the finite thermal conductivity has a 
strong destabilizing effect. Terquem & Papaloizou (1996) 
and Papaloizou & Terquem (1997) considered the stability 
of an accretion disk with an embedded toroidal magnetic 
field. These authors find that the disks containing a purely 
toroidal field are always unstable and obtained spectra of 
unstable modes in the local approximation. They argue 
that one type of modes is driven primarily by buoyancy, 
while the other is driven by shear independently of the 
magnetic configuration. 

Stability properties of purely poloidal magnetic con- 
figurations have also been well-studied. Since the papers 
by Wright (1973) and Markey & Tayler (1973, 1974), it is 
known that a poloidal field is subject to dynamical insta- 



2 



A.Bonanno & V.Urpin: Stability of the magnetic field 



bilitics in the neighbourhood of neutral points/lines if the 
field lines are closed inside the star. These authors recog- 
nizes that the magnetic field in the neighbourhood of a 
neutral line resembles that of a toroidal pinched discharge 
that is known to be unstable. Although instabilities in- 
volving significant displacements in the direction of grav- 
ity were strongy inhibited, other instabilities were not af- 
fected. The instability of poloidal configurations is rather 
fast: its growth time can reach a few Alfven crossing time 
(Van Assche et al. 1982; Braithwaite & Spruit 2006). With 
numerical simulations, the stability of poloidal magnetic 
configurations has been studied by Braithwaite & Spruit 
(2006), who apply the results to the internal magnetic con- 
figuration of neutron stars. Note, however, that a toroidal 
field might exert a stabilising influence on the instabilities 
of a poloidal field in the neighbourhood of neutral points 
(Tayler 1980). 

On the contrary, the addition of even a relatively weak 
poloidal field alters the stability of the toroidal field sub- 
stantially. For example, as first shown by Howard & Gupta 
(1962; see also Knobloch 1992; Dubrulle & Knobloch 
1993), a necessary (but not sufficient) condition for the 
instability of a toroidal field in the presence of the axial 
field reads 



> 0. 



Howard & Gupta (1962) argue that, for a fixed value of 
to, the growth rate of instability caused by condition (1) 
must vanish in the limit of a vanishing axial magnetic 
field, thereby providing a connection with the stability 
criterion obtained by Tayler (1973). Note that the pres- 
ence of a radial field is also crucial for stability prop- 
erties of rotating magnetic configurations (Bonanno & 
Urpin 2006). It turns out that configurations contain- 
ing both toroidal and poloidal fields are more stable than 
purely toroidal or purely poloidal ones (Prendergast 1956; 
Tayler 1980). With numerical simulations, Braithwaite & 
Nordlund (2006) studied the stability of a random initial 
field in the stellar radiative zone. The star was modeled 
on a Cartesian grid, and the authors found that the sta- 
ble magnetic configurations generally have the form of tori 
with comparable poloidal and toroidal field strengths. 

In the present paper, we address the stability prop- 
erties of magnetic configurations by considering the sta- 
bility of the toroidal magnetic field with respect to ax- 
isymmetric perturbations in the presence of axial fields 
of a various strength. We show that the instability may 
occur in such magnetic field configurations under the con- 
ditions that differ substantially from those imposed by the 
Tayler criterion or the necessary condition (JTJ) . We argue 
that, in some cases, the instability is caused by the new 
type of MHD waves with the growth rate proportional to 
y/B z B v where B z is the axial magnetic field. Depending 
on the profile B v (s) and the ratio B z /B v , the instability 
can occur in two regimes that have substantially differ- 
ent growth rates. We also show that the range of unstable 



wavelength in the z-direction can be essentially different, 
depending on the B z /B v ratio. 

2. Basic equations 

Let us consider the stability of an axisymmetric cylindrical 
magnetic configuration assuming plasma to be conducting 
perfectly. We work in cylindrical coordinates (s, ip, z) with 
the unit vectors (e s , e v , e z ). The inner and outer radii 
of the configuration are s\ and S2, respectively (si can 
be equal to 0). We assume that the axial magnetic field 
B z and the azimuthal field B v depend on the cylindrical 
radius alone: B z = B z (s) and B v — B v (s). Usually, the 
axial field in our model can even change the sign at some 
point so within (si, S2), S2 > sq > s\. Such a dependence 
is considered in order to mimic the stability properties of 
stars where the axial field component can generally change 
the sign. 

The equations of incompressible MHD are 

^L + (vV)v = -^ + -L(V x B) x B, 
at p Anp 



V ■ v = 0, 
3B 

— Vx(dxB) = 0, 

dt y ' ' 



(1) V-B = 0. 



(2) 
(3) 
(4) 
(5) 



It is assumed that gas is in hydrostatic equilibrium in the 
basic state and that there is no rotation, so 



Yl = J_ (V x B) x B. 

p Airp 



(6) 



In this paper we consider the stability of axisymmetric 
perturbations. Since the basic state is stationary, the time 
dependence of perturbations can be taken in the form 
exp at. Small perturbations will be indicated by subscript 
1, while unperturbed quantities will have no subscript. 
The linearized MHD equations are 

av x = -— + 7^-[(V x Bi) x B + (V x B) x BJ. (7) 
p iirp 



crBi - (B ■ V)«i + (»i • V)B = 0, 
V • Bi = 0. 



(8) 
(9) 
(10) 



Equations. (7)-(10) are homogeneous in z and admit solu- 
tions in the form of waves in the z-direction, cx exp(— ikz) 
where k is the wavevector. In stellar conditions, this anal- 
ysis applies if k satisfies the condition kH > 1 where H 
is the lengthscale of the basic state in the z-direction. 
Eliminating all variables in favour of the radial velocity 
perturbation v± s , we obtain the following differential equa- 
tion 



(a 2 +u 2 A ) 



d ( dv\ 



ds \ ds 
2fc 2 w|(l - 



a )vu 



-k v ls 



--Sujj 



4k 2 aj\ujg 



{a 2 
ds 



Vis 



^ ,(11) 
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where 



91ns 



9 In B z 
d In s 



.(12) 



In the case 5 = (B z =const), Eq. (11) recovers the 
equation derived by Acheson (1973) and Knobloch (1992). 
Once the radial velocity is known, one can calculate the 
perturbations of other quantities. For example, a pertur- 
bation of the vertical field can be expressed in terms of 
wis as 



Biz = -(B z /as)d(sv ls )/ds - v u 



dB z 
8s 



With appropriate boundary conditions, Eq. (11) allows 
the eigenvalues a to be determined. In this study, we 
choose the simplest boundary conditions and assume that 
the radial velocity vanishes at s\ and s 2 . 

3. Analytical consideration of Eq. (11) 

In this section, we consider few particular cases of Eq. (11) 
when it allows an analytical solution. The analytical con- 
sideration of simple cases can be a useful guide in un- 
derstanding the stability properties of more complex field 
configurations as well as in the interpretation of numerical 
results. 

3.1. Instability in the case B^p oc s. 

Let us initially consider the case of the azimuthal field 
being proportional to the cylindrical radius, B v = 
Bipo{s/si), where B v o is the field strength at s = s\. 
Then, the quantities B v /s and u>b are constant, and the 
last term on the left hand side of Eq. (11) vanishes. For 
the sake of simplicity, we also assume that the axial field is 
constant. After some algebra, Eq. (11) can be transformed 
into 

^ ' 1 ^ +d-i)=0, (13) 



+ 



where 
£, = ksVF 



F = 



4q^| 

{° 2 + u\y 



- 1. 



(14) 



The solution of this equation can be represented in terms 
of the Bessel functions of the order of 1 (see, e.g., Morse 
& Feshbach 1953), 



(15) 



with C\ and C 2 being constant. For the chosen boundary 
conditions, the eigenvalues of Eq. (13) are determined by 
the equation 

- n(£i)-MKi) = °> (16) 

where p = s 2 /si and £1 = kai^/F. The roots of Eq. (16) 
are real and simple. If p > 1, the asymptotic expansion of 
the nth zero is 

(„) nnsi 3As 
Si ~ ~ r - r ••• 



As 8wrs 2 



As 



si - si 



(17) 



(see, e.g., Olver 1970). Since £^ = ksiVF, we obtain the 
dispersion relation for a in the form 

(a 2 +uo\f =4/xa£w| ) (18) 

where /j, = k 2 /(k 2 + k 2 s ) and k s = Ci"V s i- Then ; 

a 2 = -ijJ 2 A ± (19) 

One of the roots is positive if the toroidal field satisfies 
the inequality 



b,m > b. (H) 



(20) 



where B v (si) is the strength of the toroidal field at the 
inner boundary. If condition (20) is fulfilled, then two ape- 
riodic modes exist with 



cri i2 w ±h 1/4 \/2uj a u b oc y/B z B v . 



(21) 



One of these modes is always unstable, and another one 
is stable and decays exactly on the timescale on which 
the unstable mode grows. The other two modes that exist 
in the considered magnetic configuration are also stable 
but oscillatory. Under condition (22), their frequencies are 
given approximately by 



cr 3i4 Ri ±i/j, 1/4 yj2cu A uJB, 



(22) 



and are also oc ^jB~ z ~B~ p . To the best of our knowledge, 
the type of MHD wave with dispersion equations (21) and 
(22) has not been considered in the literature yet. These 
modes exist in fluid only if both the axial and toroidal 
magnetic field are non- vanishing since their frequency (or 
growth/decay timescale) is proportional to ^jB~ z B~ p . Note 
that if the toroidal field is weak and condition (20) is not 
satisfied, the modes given by Eq. (18) transform into or- 
dinary Alfven waves. 

3.2. Instability of short wavelength perturbations 

If the wavelength of perturbations is shorter than the char- 
acteristic lengthscale of unperturbed quantities, then one 
can neglect the terms of the order 1/s compared to d/ds 
in Eq. (11). In this case, Eq. (11) can be transformed into 



+ 



25u)\ 



dv ls 



ds 2 s(cr 2 +u A ) ds 



k 



1 



M<4 



+ 



2w|(l 



(a 2 + , 
a) 



2 ^,2 
A) 



{a 2 



Vis = 0. 



(23) 



In a short-wavelength approximation, the coefficients of 
this equation can be treated as constant, and the solution 
can be taken in the form F oc exp(— ik s s) where k s is the 
wavevector in the radial direction. The dispersion relation 
corresponding to Eq. (23) is 

(a 2 + u A )(a 2 +uj\ + 2A) - ^u; 2 a lj 2 b = 0, (24) 



where 



(1 — a)fj,oj 



ik 
sq- 



2 SoJ A , 



// 



Q 2 



Q 2 = k 2 + k 2 
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The solution of Eq. (24) is 



2 2 

a =-u A 



-A±JA 2 +4/iu) 



2/ i 2 
A B ' 



(25) 



If the azimuthal field is vanishing (wg =0), then the roots 
are 



'1,2 



= — UJ 



A . 



'3,4 



2z 

^Q 5 



The first two modes are always stable. The modes 3 and 4 
are usually stable since the second term in the expression 
for (T 3; 4 is of the order of l/(sk s ) <§C 1 compared to the first 
term and should be neglected in a short-wavelength ap- 
proximation. The only exception is the region in a neigh- 
bourhood of the neutral point (or line) where B z (s) — > 
and, hence, uja — > 0. For example, if -B z (s) goes to zero 
oc (s — so), then the first term in erf 4 is oc (s — so) 2 , but 
the second one is oc (s — s ). Therefore, the region always 
exists in the neighbourhood of so where the second term 
dominates the first one. In this region, the roots 3 and 4 
are approximately given by 



03,4 




OJA- 



(27) 



One of the roots (27) corresponds to the unstable mode. 
Therefore, poloidal magnetic configurations with the neu- 
tral point (line) are always unstable in the neighbourhood 
of this point (line). This fact was first pointed out by 
Tayler (1973). 

However, the presence of the azimuthal field can drasti- 
cally change the stability properties of the magnetic field 
even in a neighbourhood of the neutral point. Near the 
neutral point, the azimuthal field is stronger than the ax- 
ial one (if the former does not have zero at the same point 
s ), then from Eq. (25) we have with the accuracy in terms 
of the order of uj a 



a\ 2 ~ 2(a — 1)/mj 



2 

B ! 



; 3,4 



1 



a 



(28) 



It turns out that the stability is determined by the prop- 
erties of the azimuthal field rather than the presence of 
a neutral point in the magnetic configuration. If a > 1, 
then one of the modes (1, 2) is unstable. On the contrary, 
one of the modes (3, 4) is unstable if a < 1. Note that 
the instability of the modes 3 or 4 is less efficient in the 
neighbourhood of the neutral point than in the rest vol- 
ume (see below) because uja is small there and, hence, the 
growth time of instability is long. 

Consider now the local stability of a region that is far 
off the neutral point. In this case, the term in A propor- 
tional to S can be neglected in Eq. (25). Then, we have 



a + a 2 cr + a = 0, 



where 



(29) 



a 2 = 2uj 2 a + 2[MjJ%(1 - a), a = uj a [lu a — 2pno\(\ + a) 



Equation (29) has positive (unstable) roots if either 02 < 
or a < 0. These two conditions are equivalent to 



a > 1 + (ksx)' 



and 



(26) a>-l + (fcs 1 ) ; 



Eft 



(30) 



(31) 



respectively. Obviously, if condition (31) is fulfilled, then 
condition (30) is also satisfied. Therefore, the true crite- 
rion of instability is given by Eq. (31). Depending on the 
parameters, the sufficient condition of instability (31) can 
differ substantially from the necessary condition (1) that 
is equivalent to a > — 1. The roots of Eq. (29) are 



■■-uj a - 



M <4(1 - a) ± ^^4(1-^)2+4^1(32) 



In the case of a weak axial field (ub > uja or B v > 
(ks\)B z ), the behaviour of instability is rather different 
for a > 1 and a < 1. If a > 1, the unstable root of 
Eq. (32) is equal to 



a w Ld B y/2Li(a - 1). 



(33) 



In this case, the instability grows on the timescale of the 
order of the time it takes an Alfven wave in the field B v 
to travel to the distance s\. This regime recalls the insta- 
bility of a purely toroidal field that occurs only at a > 1 
and grows on the same timescale ~ (Braithwaite & 

Nordlund 2006). If a < 1, then the root corresponding to 
instability is given by 



1 



-OJA- 



(34) 



For such profiles of the toroidal field, the instability essen- 
tially grows slower. Its growth time is of the order of an in- 
verse Alfven frequency associated with the axial magnetic 
field, B z . This instability is less efficient in a relative neigh- 
bourhood of the neutral point where uj a is small. Note 
that the value a w 1 distinguishes between two regimes 
only in the limit of weak B z . If B z is stronger, then the 
distinguishing value becomes higher. 

Note that, in the case a = 1, Eq. (32) recovers the 
dispersion relation (19). If the axial magnetic field is suf- 
ficiently strong, axisymmetric perturbations are stable for 
any profile B v {s). 

4. Numerical results 

In this section we discuss the stability of the magnetic 
configurations containing axial and azimuthal fields by nu- 
merically solving the eigenvalue problem. In calculations, 
we assume that the dependence of the azimuthal magnetic 
field on s is given by 



B lr 



B v a — 
1 si 



(35) 
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where B v o is the field strength at the inner boundary. To 
calculate the growth rate of the instability, it is convenient 
to introduce dimcnsionless quantities 



x=—, q = ks!, T= , u B o = — 

si ubo si 

Then, Eq. (11) transforms into 



^ (36) 



B, 



d f dvi 



dx \ dx 
2(1 -a) 



X 



T 2 + q 2 e 2 



<r < 1 - x 2 ^' 

25 



2J1 



(r 2 + g 2 e 2 ) 2 



Vis = - 



q 2 e 2 



x T 2 + q 2 e 2 \ dx 



+ 



• (37) 



The parameter e can depend on s in this equation since 
B z is generally a function of s in our model. The in- 
ner and outer boundaries correspond to x\ = 1 and 
x 2 = 2, respectively, where we have vi s (xi) = vi s (x2) = 0. 
Equation (36), together with the given boundary condi- 
tions at the extrema, is a two-point boundary value prob- 
lem that can be solved by using the "shooting" method 
(Press et al. 1992). In particular in order to solve Eq. (36), 
we used a fifth-order Runge-Kutta integrator embedded 
in a globally convergent Newton-Rawson iterator. As the 
spectrum is discrete, we first found the lowest eigenvalue 
for the analytical solution (19), and then we gradually 
changed q, B z /B v and a in order to explore the parame- 
ter space. In this way we checked that the eigenvalue was 
always the fundamental one, as the corresponding eigen- 
function had no zero except the one at the boundaries. 
Higher radial wavenumbers corresponds to lower growth 
rates, as it is also apparent from Eqs.(20)-(21). 

The radial dependence of the eigenfunctions is shown 
in Fig. 1 for a = 2 and B z /B v = 0.5 {5 = 0). The curve 



corresponding to q = 1 describes the unstable mode, but 
the other 3 curves represent stable modes. The eigenfunc- 
tions for q < 1, representing unstable modes as well, al- 
most coincide with the curve q — 1 and are not shown 
in Fig. 1. Note that, for other values of the parameters, 
the unstable eigenfunctions can also have sharp maximum 
near the outer (if a > 0) or inner (if a < 0) boundaries. 
It seems that the maxima tend usually to be located near 
the boundary with the strongest field. 

In Fig. 2, we plot T = u/lobo 

function of q = ks\ 

for a weak constant axial field. The addition of even very 
weak B z changes the stability properties substantially. A 
purely toroidal field should be unstable if a > 1. Our 
calculations show, however, that even the profile of the 
toroidal field with a = is unstable in this case de- 
spite a very low energy of the axial field compared to 
that of the toroidal field (<~ 1%). It is also worth noticing 
an important difference between the necessary condition 
of instability (1) and the sufficient condition. In accor- 
dance with inequality (1), the necessary condition requires 
a > a c = — 1, but our calculations show that the insta- 
bility occurs at essentially larger a c ~ —0.1. The growth 
rate as a function of q always has a smoothed maximum at 
q ~ 5 — 20 depending on a. The value of q corresponding 
to the maximum increases with increasing a. The figure 
clearly illustrates the difference between the two regimes 
of the instability discussed in Sect. 4. If a < 1, the insta- 
bility is rather weak and grows on the AlfVen timescale 
characterised by the axial field. In contrast, if a > 1, the 
instability is much more efficient, and the growth time is 
close to the Alfven timescale for the toroidal field. 
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The dependence of T 2 on q for a stronger axial field 
is shown in Fig. 3. An increase of the axial field makes 
the magnetic configuration more stable. For example, in 
the case a = 2, the maximum growth rate of instability is 
w 2.5 uobo if e = 0.1, but it only becomes w 1.7 wso if £ = 
0.5. As mentioned, the critical value of a that determines 
the onset of instability is ~ —0.1 in the case e = 0.1. If 
£ = 0.5, then instability can occur only if the toroidal field 
increases enough rapidly with s, and the critical value a c is 
close to 1. In a stronger axial field, the transition between 
two regimes of instability (see Eqs. (30) and (31)) occurs 
at a higher value of a as was predicted in Section 3.2. 

In Fig. 4, the normalized growth rate is shown for 
B z /B v = 1. Such a strong axial field stabilises the mag- 
netic configuration drastically. The instability occurs only 
if the toroidal field increases very rapidly with s (a > 2). A 
qualitative behaviour of unstable eigenvalues is the same 
as in Fig. 2 and 3, but the maximum growth rate corre- 
sponds to lower values of q. 

In Fig. 5, we plot the dependence of T 2 on £ = B z /B v0 
for q = 1. The profiles of the toroidal field with a = 0.5 
and 1 are stable if B z = 0. Therefore, these curves start 
with r 2 = 0. In contrast, the toroidal field profile with 
a = 2 is unstable even if the axial field is vanishing and, 
therefore, T 2 (£ = 0) ^ 0. All the curves shown in Fig. 5 
exhibit a qualitatively similar behaviour: the growth rate 
reaches a rather flat maximum at B z /B v0 ~ 0.3 — 0.5 
and then decreases fast enough. The quantity T 2 becomes 
negative and, hence, the magnetic configuration is stable 
if B z is sufficiently strong. The higher is the value of a, 
the stronger the axial field required for stabilization. 




0.0 0.5 1.0 1.5 

B /B 

Fig. 5. The dependence of T 2 on B z /B v v = 1 for q = 1 
and different values of a. 

In Fig. 6, we plot the dependence of T 2 on q for the 
axial magnetic field that depends on s. We choose this 
dependence in the form 

B z (s)=B z0 (l--)=B z0 (l-—), (38) 
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where xq — sq/si, s is a radius at which B z changes sign, 
S2 > so > si. We have for such B z 



x 



x — x 



Xq — X 

x 



(39) 



where eo = B z o/B v q. Calculating models in Fig. 6, we 
suppose x = 1.75 and e = 0.6. The ratio B z /B v o aver- 
aged over s is equal to 0.1 in this case, thus the results can 
be compared with those in Fig. 2 where B z =const and 
B z /B v q = 0.1. Although the axial field has the neutral 
line at so = 1-75, the stability properties are not very dif- 
ferent from the similar configuration with no neutral line 
(Fig. 2). In Fig. 6, the critical value of a that distinguishes 
between stable and unstable configurations is close to 0. If 
the axial field has no neutral line (Fig. 2), then the critical 
value is fairly close to this value (~ —0.1). The maximum 
growth rate is of the same order of magnitude for the same 
a. The only difference is that configurations with the neu- 
tral line (Fig. 6) are unstable for a wider range of q. 

In Fig. 7, we plot T 2 versus q for the case of a stronger 
axial field (B^/B^o = 3) with the same neutral line at 
s = 1.75. The quantitative behaviour is similar to what is 
depicted in Fig. 6, although the growth rates are lower be- 
cause of the higher value of B z0 /B v0 - The averaged ratio 
B z /B v0 is 0.5 in this case, and the results can be com- 
pared to Fig. 3 where B z /B v q = 0. It appears that the 
presence of a neutral line does not change the stability 
properties substantially. For instance, the critical value of 
a is ~ 1 for both configurations. The maximum growth 
rate at a given a is of the same order of magnitude as well. 
However, the rage of unstable q is wider for configurations 
with the neutral line. 



5. Discussion 

We have considered the hydromagnetic stability of cylin- 
drical configurations containing the toroidal and axial 
magnetic fields. Dissipative effects were neglected in our 
study. We treated a linear stability assuming that the be- 
haviour of small perturbations is governed by equations 
of incompressible hydrodynamics. This approximation is 
justified if the magnetic field is subthermal and the Alfvcn 
velocity is low compared to the sound speed. The stability 
of the magnetic configurations is a key issue for under- 
standing the properties of various astrophysical bodies, 
such as peculiar A and B stars, magnetic white dwarfs, 
neutron stars, etc. Even though various dynamo models 
predict that the toroidal field should be typically stronger 
than the poloidal one, the effect of a poloidal field on the 
stability usually cannot be neglected. 

To demonstrate this, we treated the simplest model of 
a highly conducting fluid between two cylindrical surfaces. 
We assumed that the toroidal and axial fields depends 
on the cylindrical radius alone. In a short-wavelength ap- 
proximation, we derived the growth rate and a sufficient 
criterion of instability analytically (Eqs. (30)-(31)). For 
large-scale perturbations, the condition of instability and 
its growth rate were calculated numerically. The analytical 
and numerical results are in good qualitative agreement. 
The obtained conditions of instability differ substantially 
from what is predicted by the necessary condition (1). For 
instance, according to Eq. (1), if the instability occurs in 
the magnetic configuration, then the toroidal field profile 
satisfies the condition a > — 1. In fact, the instability oc- 
curs only if the toroidal field decreases with s much slower 
(or even increases): the critical value of a is w —0.1 if 
B z /B v0 = 0.1 and w 1 if B z /B ip0 = 0.5. If B z depends 
on s, then the critical values of a should be even higher. 
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Depending on the profile of the toroidal field and the 
strength of the axial field, the instability can arise in two 
essentially different regimes. In the case of a weak axial 
field, B v o ^> B z , the value of a that distinguishes between 
the regimes is ss 1. If a > 1, then the instability grows 
on the Alfven timescale determined by the toroidal field 
and is rather fast. If a < 1, then the instability is much 
slower and grows on the timescale determined by the axial 
field. The transition between these two regimes occurs at 
larger a if the axial field increases. The efficiency of the 
considered instability turns out to be rather low if a < 1 
and B z is weak. 

It is worth noticing the very particular properties of 
instability in the case a ~ 1. In such a configuration, 
the particular type of MHD waves is given by the disper- 
sion equation (19). The growth rate of these waves (or 
the frequency, if the waves are stable) is proportional to 
the product of B z and B v . These waves cannot exist in 
purely toroidal or purely poloidal fields. The instability of 
the configuration with a » 1 is caused by the generation 
of this particular type of wave. These waves can probably 
determine the instability of magnetic configurations near 
the axis of symmetry where B v — > 0. 

A sufficiently strong axial field always suppresses the 
instability. For more or less plausible values of a < 1, the 
strength of the axial field stabilising the configuration is 
~ 0.1 — 1B V . A much stronger field is required, however, 
to stabilise the configuration with larger a. 
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